home *** CD-ROM | disk | FTP | other *** search
/ NetNews Offline 2 / NetNews Offline Volume 2.iso / news / comp / lang / c-part2 / 12349 < prev    next >
Encoding:
Internet Message Format  |  1996-08-05  |  2.1 KB
  1. Path: news.bu.edu!usenet
  2. From: Xiangyang Wang <xianwang@chem.bu.edu>
  3. Newsgroups: comp.lang.fortran,comp.lang.c
  4. Subject: Re: Efficient Column-wise FFT Routines?
  5. Date: Sun, 31 Mar 1996 12:19:52 -0500
  6. Organization: Boston University
  7. Message-ID: <315EBEB8.41C67EA6@chem.bu.edu>
  8. References: <warrenp-2603961110500001@ren.colorado.edu>
  9. NNTP-Posting-Host: radon.bu.edu
  10. Mime-Version: 1.0
  11. Content-Type: text/plain; charset=us-ascii
  12. Content-Transfer-Encoding: 7bit
  13. X-Mailer: Mozilla 2.0 (X11; I; SunOS 4.1.4 sun4m)
  14.  
  15. Pete Warren wrote:
  16. > 
  17. > Greetings all,
  18. > 
  19. > I am trying to perform large scale FFT's on complex data in a column (or
  20. > row) wise formats.  Does anyone know of a good, proven routine for either
  21. > Fortran (preferred) or C ?  There are routines for complex 1-D data and
  22. > for real column wise data but not for both.
  23. > 
  24. > What I'm actually trying to do is build a code interface node for Labview
  25. > that performs this function efficiently.
  26. > 
  27. > I would appreciate any help or pointers.
  28. > 
  29. > Thank you,
  30. > 
  31. > - Pete Warren
  32. > warrenp@colorado.edu
  33. The following fft subroutine is from web resourse( I am not sure where)
  34. , maybe it helps.
  35.       SUBROUTINE FOUREA( N, ISI)
  36. C
  37. C THE COOLEY-TUKEY FAST FOURIER TRANSFORM IN ANSI FORTRAN
  38. C
  39. C DATA IS A ONE-DIMENSIONAL COMPLEX ARRAY WHOSE LENGTH, N, IS A
  40. C POWER OF TWO.  ISI IS +1 FOR AN INVERSE TRANSFORM AND -1 FOR A
  41. C FORWARD TRANSFORM.  TRANSFORM VALUES ARE RETURNED IN THE INPUT
  42. C ARRAY, REPLACING THE INPUT.
  43. C TRANSFORM(J)=SUM(DATA(I)*W**((I-1)*(J-1))), WHERE I AND J RUN
  44. C FROM 1 TO N AND W = EXP (ISI*2*PI*SQRT(-1)/N).  PROGRAM ALSO
  45. C COMPUTES INVERSE TRANSFORM, FOR WHICH THE DEFINING EXPRESSION
  46. C IS INVTR(J)=(1/N)*SUM(DATA(I)*W**((I-1)*(J-1))).
  47. C RUNNING TIME IS PROPORTIONAL TO N*LOG2(N), RATHER THAN TO THE
  48. C CLASSICAL N**2.
  49. C AFTER PROGRAM BY BRENNER, JUNE 1967. THIS IS A VERY SHORT VERSION
  50. C OF THE FFT AND IS INTENDED MAINLY FOR DEMONSTRATION. PROGRAMS
  51. C ARE AVAILABLE IN THIS COLLECTION WHICH RUN FASTER AND ARE NOT
  52. C RESTRICTED TO POWERS OF 2 OR TO ONE-DIMENSIONAL ARRAYS.
  53. C SEE -- IEEE TRANS AUDIO (JUNE 1967), SPECIAL ISSUE ON FFT.
  54. C
  55. -Xiangyang
  56. .
  57.